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ABSTRACT 

The  sensitivity  of  planetary  waves  to  various  initial 
conditions  and  thermal  forcing  is  examined  using  a  linear 
global  primitive  equation  spectral  model. 

Initial  conditions  are  ottained  by  switching  on  a 
analytic  heat  source  and  then  integrating  the  model  equa- 
tions out  to  30  days.  An  averaging  procedure  is  used  to 
eliminate  any  transient  modes  vhich  remained  after  integra- 
tion so  that  the  initial  conditions  represent  steady  state 
solutions.  Additional  integrations  are  performed  in  which 
'errors'  are  introduced  into  the  forcing  and  initial 
conditions. 

Results  of  the  study  are  examined  using  polar  phase  vs 
amplitude  plots  (harmonic  dials)  of  various  spherical 
harmonics.  Results  indicate  that  planetary  waves  are  not 
sensitive  to  errors  in  the  forcing  or  initial  conditions. 
However,  this  lack  of  sensitivity  is  most  likely  due  to  the 
simplified  initial  conditions  used  in  the  model  rather  than 
to  the  inherent  sensitivity  of  planetary  waves. 
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I.  INT1QPJDCTI0N 

Atmospheric  predictability  studies  (Lorenz,  1969)  indi- 
cate each  scale  of  motion  has  a  limit  to  it's  predict- 
ability. Small  scale  motions  are  theoretically  predictable 
to  an  hour,  synoptic  scales  to  a  few  days  and  planetary 
scales  to  a  few  weeks.  However  a  number  of  studies,  Lambert 
and  Merilees  (1978),  Baumhefner  and  Downey  (1978),  and  Morse 
(1983),  have  shown  that  synoptic  scale  motions,  rather  than 
planetary  scale  motions,  are  the  most  accurately  forecast. 
Fig.  1.1  is  taken  from  Daley  et  al  (1981).  It  gives  the 
500mb  geopotential  forecast  error  from  the  6- level  National 
Center  for  Atmospheric  Research  global  circulation  model  as 
a  function  of  wavenumter.  The  error  for  each  wavenumber  has 
been  normalized  by  the  natural  atmospheric  variance  for  that 
wavenumber.  Note  that  as  previously  stated  the  synoptic 
scales  (wavenumbers  4-8)  are  in  fact  the  most  accurately 
forecast.  This  result  does  not  seem  to  be  dependent  on  the 
type  cf  model  used.  Baumhefner  and  Downey  examined  a  number 
of  different  models  and  found  similar  results.  The  error  in 
forecasting  planetary  scale  waves  might  not  seem  to  be  that 
important  given  that  the  error  in  the  synoptic  scales  is 
smaller  and  that  these  are  the  scales  of  motion  which 
produce  most  of  the  day  to  day  weather  changes.  The  impor- 
tance of  an  accurate  planetary  scale  comes  to  light  when  one 
starts  to  consider  medium  to  long  range  forecasts  (up  to  1 0 
days).  Since  the  planetary  waves  often  act  to  steer  the 
smaller  synoptic  scale  disturbances  an  improved  planetary 
scale  forecast  would  presumably  lead  to  an  improved  forecast 
on  the  synoptic  scale.  In  addition,  the  planetary  waves 
contain  a  major  portion  of  the  eddy  kinetic  energy  in  the 
atmosphere  and  again  a  better  long  wave  forecast  would  lead 
to  a  better  synoptic  scale  forecast. 


There  are   at  least  two   possible  reasons   why  planetary 
waves  are  not  forecast  as  well   as  theory  suggests.    One  is 
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Figure  1-1   Forecast  Errors  for  1r  2 
and  3  Days  of  Integration. 


that  the  model  dynamics  for  planetary  waves  are  inadequate. 
Given  that  planetary  waves  are  quasi-stationary  in  nature 
and  are  forced  in  part  by  differential  heating  and  orography 
it  would  seem  reasonable  to  suspect  that  inadequacies  exist 
in  the  dynamics  of  the  model's  forcing.  Another  possibility 
is  that  current  initialization  procedures  are  not  adequate 
for  planetary  scales. 

It  is  the  purpose  of  this  study  to  examine  the  nature  of 
planetary  waves   and   to   deternine   the   sensitivity   of 
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planetary  waves  to  initial  conditions  and  forcing.  Chapter 
two  of  this  thesis  describes  the  model  used  in  this  study. 
The  third  chapter  outlines  the  nature  of  the  experiment.  In 
the  fourth  chapter  the  results  of  the  experiment  are 
presented  and  in  the  final  section  these  results  are 
discussed  and  some  conclusions  are  drawn. 


1  1 


II.  MODEL  DESCRIPTION 

The  model  used  in  this  study  is  a  baroclinic  spectral 
transform  model  which  was  developed  by  Rosmond  (1977). 
Similar  formulations  have  been  carried  out  by  Hoskins  and 
Simmons  (1975)  and  by  Bourke  (1974).  The  model  is  config- 
ured to  include  friction  and/or  diabatic  heating.  The 
specifics  of  how  friction  and  diabatic  heating  are  included 
in  the  model  will  be  discussed  in  the  next  section.  The 
basic  equations  of  the  model,  in  sigma  coordinates,  are  as 
follows: 


||  =  -V.(C  +  f)V-k^x(RTVVa^)+k.Vx?  <2'1J 

3 1  ^       do 


~   =  k-Vx(;  +  f)V-V-(RTVq  +  a^-)-V2(4)+^  ) +7  •  F 

d  t  n    dO  2 


(2.2) 


||  =  -D-V.Vq-ii  (2.3, 

il    _£1  (2.5) 

da    '      a 

where : 

C  -  vorticity 

D  -  divergence 

1  -    temperature 

8  -  potential  temperature 

tt  -  surface  pressure 

V  -  horizontal  velocity  vector 

<t>  -  geopotential  height 

R  -  gas  constant 

Cp  -  specific  heat  at  constant  pressure 

f  -  Coriolis  parameter 
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a  -    vertical  coordinate  (a  =  P/tt) 

o  -    vertical  velocity    (6  =  -~  ) 

g   -    In 

P  -    p 

^    -    R/Cp 

F  -  frictional  force 

The  continuity  equation  (Eg.  2.3)  may  be  rewritten  by 
integrating  with  respect  to  signa  and  applying  the  boundary 
conditions  d(0)=a(1)=0.  Thus  the  integral  of  Eg.  2.3  maybe 
written  as 

|f   =    -D+G  (2-6) 

where    (      )    =  fH      )da  and    G=    V  -Vg. 

The  vertical  velocity,  j,  may  be  obtained  diagnostically 
by  substituting  Eg.  2.6  into  Eg.  2.3  and  integrating  in  the 
vertical  to   obtain 


d    =    (D+G)o-/°(D+G)da  (2"?) 

which   uses      5  (0) =0 

The   first    law    of    thermodyna nics    Eg.    2.4   can    be    written 

|1  =  4-VT-oKd^    (Ta-K)+KT(|a+V.?q)4  (2-8> 

dt  9a  3t  CJp 

In    order   to  apply    semi-implicit   differencing    it    is    necessary 
to    divide   temperature  as    follows: 


T  =  T*(a)+T!(a,  X,  Y,  t)  (2"9) 

where  T*  represents  an  appropriately  averaged  temperature. 

The   basic  equations   can  be   conveniently  written  in 
spherical  coordinates  by  defining  the  following  operator: 
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a(a,b)   =       *     |a*gb  (2.10) 

1-x*    9A   9x 

Using    Eqs.       2.9    and    2.10    the      basic    equations   can    be    written 
as    follows: 

||=-a(A,B)  (2-11) 

ut 

||  =   -a(A,B)-?2(E+$+RT*q)  (2*  12) 

U  =  -a(UT,,VT,)+DT,-a<a^-<Ta"lc)+icT(G-G-D)  {2'  13) 

dt  9a 

|f--(D+G)  f2-11" 


where 

A  =    (;+f)U+a|X+C^)(l-x^)|^2FT 

1-x2   9X     v9x 


E  =   (U2+V2)/2(l-x2) 
U  =  uCcosWr 

V  =  vCcosWr 
x  =  sinY 

V  -   latitude 
X  =  longitude 


HI 


Egs.  2.10  -  2.15  are  the  basic  equations  used  in  the  model. 
These  eguations  are  represented  spectrally  in  the  horizontal 
and    finite    differenced    in    the    vertical. 

A.       VERTICAL    STRUCTOBE 

The  vertical  structure  of  the  model  follows  the  develop- 
ment given  by  Arakawa  and  Suarez  (1983).  The  variables  are 
staggered  in  a  so  that  X, ,  D,  U,  V  and  T  are  carried  at  the 
mid-point  of  each  layer  where  a  =  a^  and  a  is  carried  at  the 
top    and     bottom    of    each    layer      where  o    =  a,     .  The    vertical 

structure  is  illustrated  in  Fig-  2.1.  The  finite  difference 
form    of    Egs.    2.10   -    2.15    are 


If  =  -«k<A,B) 


l^k    =    a,  (B^A)-^2^  +RT*q+En  ) 


3t 

i£    - 

8t    " 


'k       ' 
-(G+D) 


k 


r-k    =    -a.  (UT*    VT'  )T'D1  +    T*(G1-G-D)- 


3t 


'k 


k         k      k 


^k+lBJ<(^+Tk+1-Tk>+    W   f^T^+Q/Cp 

v*k+1  =cp(pk+i-pk>^¥k+^> 

k  k+1 


LM' 


LM 
k 


av+i    =    a,  .,(G  +  D)-5       (G.+D.)Aa. 
k+1  k+1  j  =  l      3       ]         3 


where 


pk  =  (ak+1*)K 


P      -      <  V  \  ( ^1  +  K      1  +  K  \   /  / 

P    -    (ITK-)(ak  +  l0k      )/(ak  +  l-ak} 


(2.  16) 
(2.17) 
(2. 18) 
(2.19) 


(2.  20) 
(2.21) 
(2.22) 
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CTr 

'-  CT, 

=0 

{, 

D, 

T, 

Aa, 

0, 

1 

1 

i 

3  K-1            K-1           K-1            K-1 

<Vl 

I                                                       ?k    DK    TK    AaK    0K 

X"V                              W                              -"*v 

y~         r\         t 

=>  K+1            K  +  1           K  +  1            K  +  1 

^K+1 

^K*. 

I 

I 
I 

<Tlm     Dlm     Tlm     A<tl» 

0LU 

LM  +  1            LM  +  1        r  LM  +  1 

Y/////// 

/     /     /     / 

Figure  2.1    Vertical  Structure. 
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A    =     U  +  f)Uk+{dk+1(Vk+1-Vk)+dk(Vk-Vk_1)}/2Aok 

+  {RTV(l-x2)/r2}l^   -(/T^2/r)Flw 
B    =    U+f)Vk+{dk+1(Uk+1-Uk)+dk(Uk-Uk_1)}/2Aa 


k 


-(CpRT^/r)-^   +(/T^2/r)F 


k,x'3x    '^"*    /r^rA 


G     =      Vk-^TT/7T 


D    =    v*«v\ 
k 


A,     =     (P.  -P.  >/(Pv.-,-Pv) 

k  k      k  k+1      k 


k 


=    (Pk+l-Pk)/(Pk+l-p^  =  l-Ak 


Egs.  2.16  -  2.22  can  be  written  in  matrix  form,  such 
that  the  terms  on  the  right  hand  sides  contain  all  the  terras 
which  are  to  be  evaluated  explicitly  and  the  left  hand  side 
contains  those  terms  which  are  to  be  evaluated  implicitly. 
Eg.  2.20  can  be  combined  with  an  integrated  finite  differ- 
ence   form   of   Eg.      2.  15   to  obtain 

4    =    C   T+d)  (2*  23) 

where  C  is  a  sguare  matrix  and  the  other  quantities  are 
column  vectors.  The  finite  difference  form  of  the  surface 
pressure  tendency  equation (Eg.   2.14)  is 

U  =  Z(Gv  +  Dv)Aav  (2-2tO 

3t       k   k    k 

which  can  be  written  in  matrix  form  as 


|£  =  -NT(G  +  D)  (2'25) 
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T     . 
where      N      is   the      transpose      of      a    constant      column      vector 

Similarly   Eg.      2.22    can    be    written 

o    =    Z(G+D)  {2m  26) 

The    next   to   last    term  in    Eg.     2.  19    is 

pk    rp  nr    ^  iA    A  ,m         Pv    m  ^2.  27) 


k(A'^w-Vt»A.1«V  H.i» 


k-1 

For  the  purpose  of  seni-implicit  formulation,  the  tempera- 
ture is  separated  according  to  Eg.  2.9.  The  mean  part  of 
that    term    can    be   written   as 

VlVF^  +  1-T*)+aAk_l(Tg-    St-T*)     =    M(G  +  D)  <2"28> 

K+ l  k- 1 

Egs.  2. 17,  2.19,  2.23,  and  2.25  may  now  be  written: 
!f+V2(r+R?*q>=a(B,-A)-V2E  =  ^  (2.29) 

|i+NTD  =  -nt5  <2-30) 

||+g  2  =  KT  I2-3" 

J'  =  C  f  ^2*  32^ 

T 
where  £      =    M   *      kT*N      and   <£=<}>-*       and    KD    and      KT  represent 

terms    which   have    not    been    explicitly    separated    out. 

The  semi-implicit  time  differencing  is  achieved  by  eval- 
uating the  terms  on  the  left  hard  sides  of  Egs.  2.29,  2.30 
and    2.31    implicitly.  The    remaining    terms    and      Eg.       32    are 

evaluated      explicitly      using      leapfrog      differencing.  The 

difference   eguations   can   now    be   written 

D    ^n+AtV2(C    T    .1+R^*qri4.1)    = 
n  +  1  =      n  +  1    —     ^n+1 

'     +  (2.  33) 

D      rfitV2(CT       , +RT*q       _)+2At(Kn)_ 
n-1  —      n-1    _      ^n-1  L)   n 
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fn+l+At25n+l    =    V1-Atg8n_1-2At(gT)n  (2.34) 

q.^+AtNTD    =    q       -.-AtivTD       1-2AtNiG  U*  ^' 

n+i  ^n-1  n-1  n 

Now    the    following      eguation    for    D         can    be      found    by    substi- 
tuting   Eg.    2.34    and    2.35    into    2.33: 

BA  +  1    =    B+Snil+.24t(«D)n-2At^(cfn_1+rf*qn_1+C(RT)n    <2-36' 

-RT-N    G 
=  n 

where   the   matrix   operator   B    is 


B±    =    At2(C    Q±RTN    )V2+I_  l  ' 

B.       SPECTRAL    FORMULATION 

The    eguations    (2.  10) -  (2. 15)         are    represented    spectrally 
in      the      horizontal.  The      variables      are      represented      as 

follows: 

CCA,    x,    a,    t)    =    I       .   I    ,     ,Cm(c,    t)Pm(x)eimA    =  (2.33) 

m=- J    n= j  m|     n       '  n 

J  J 

E       .   S    ,     | CmYm 
m=-J    n=    m      n    n 

m     *  m         -m 

where  C    is   some    variable   and    (Cn  )       =    (-1)  _    Cn 

m   =    zonal   wavenumber 

n   =    meridional      index,      and    n-|m|    gives    the      number    of    zeros 

between   the   poles      (-1    <   x   <    1)      of      the   associated    Legendre 

function 

J  =   truncation  limit  (for   this  study   triangular  truncation 

was  used  with  J=3) 

A-   (1-1) /2   nondimensional   zcnal    coordinate   index 

(1  <  1  <  16) 


19 


Note  that  the  separation  is  such  that  the  coefficients  C™ 
are  functions  of  time  and  the  vertical  and  the  spherical 
harmonic  Y n are  horizontal  functions  of  space.  The  normali- 
zation  and  orthogonality   properties  of   the   Ym  allow   the 

n 
coefficients  to  be  obtained  as  follows: 


Cm    =    <C,Ym>    E    -l/n2TT   /^CY^xdA  <2-39> 

n  n  Utt    0         -1         n 

The  non-linear  terms  are  computed  using  the  transform  method 
following   Haltiner      and   Williams    (1980).  The    longitudinal 

direction   is    treated      with   a    Fast    Fourier      Transform    and    the 
latitudinal    direction   uses   Gaussian    Quadrature.         The    number 
of    latitudes,    N    and    longitudes    ft   satisfy. 
N   >    3J/2    +1  M    >    3J+  1 

The  number  of  points  are  chosen  so  that  there  will  be  no 
aliasing  from      the    product    terms.  For   this  study      N=4   and 

M=16.  It      was      discovered      just      before      this      study      was 

completed  that  number  of  Gaussian  latitudes  used  to  compute 
the  latitudinal  integrals  was  actually  J+1  and  not  3J/2  +  1. 
However,  model  integrations  using  J  +1  latitudes  are  not 
qualitatively  different  than  these  using  3J/2  +  1  latitudes. 
Even  the  quantitative  differences  are  small,  and  will  occur 
mostly      in    the      non-linear    integrations.  Since    the      cases 

examined  in  this  study  are  carried  out  using  a  linear  form 
of   the    model,    one    would   not    expect    much    error. 
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III.  EXPERIMENT  DESIGN 

The  intent  of  this  study  is  to  specify  a  known  analytic 
heating  function  on  a  planetary  scale  and  then  integrate  the 
equations  until  the  model  atmosphere  reaches  a  steady  or 
quasi-steady  state.  These  steady  state  solutions  will  be 
altered  and/or  the  heating  will  be  changed.  The  equations 
will  then  be  integrated  again  to  determine  the  sensitivity 
of  planetary  waves  to  changes  in  the  forcing  and  initial 
conditions 

Heating  is  introduced  into  the  model  via  the  thermody- 
namic eguation  (2.4).  Q  is  specified  to  be  a  Legendre 
function  of  the  form: 

QCX,  x,  o)  =  l_         z    ,  ,Am(a)Pm(x)eimX  (3-1) 

m--J  n= | m|  n     n 

A   is  specified  to  have  the  vertical  structure: 
n 

n      .3 

Only  one  wave  form  is  forced  at  a  time  and  the  other  heating 
amplitudes  are  set  to  zero.  The  constant  A  is  chosen  to 
give    a      reasonable    temperature    response.  In   most      cases   a 

temperature   response    of   a   few    degrees   is    used. 

Before  the  heating  is  "tuned  on"  the  model  atmosphere 
is   specified   to   be   in  a   state    of  solid   rotation 

IT  =        cos 
V   =   0 
such    that      U   =      20    m/s(15.5      °/day)       at      the   equator.  The 

surface  pressure  field  is  chosen  to  balance  the  mean  flow. 
The  initial  vertical  temperature  distribution  is  given  by 
the   standard   atmosphere    values      consistent    with    the    pressure 
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distribution  given  by  solid  rotation.    A  lapse  rate  of  6.5° 
per   kilometer   is   specified  tc   the  top   of   the   model 
atmosphere. 


TABLE  I 
Waves  Forced  for  linear  Cases 


Case  Wave  Forced 


1 
2 

3 
4 
5 
6 


m 

n 

1 

1 

1 

2 

1 

3 

2 

2 

2 

3 

3 

3 

In  all  of  the  experiments  the  model  is  run  in  a  linear 
mode.  In  this  mode  the  amplitudes  of  all  zonal  wavenumbers 
except  those  being  forced  are  held  constant.  The  model  is 
run  in  a  linear  mode  in  these  experiments  to  allow  a  more 
definitive  interpretation  of  the  results.  Six  linear  cases 
are  integrated,  each  case  with  a  different  forced  wave 
form.  For  example  in  case  1  (see  Table  I)  zonal  wavenumber 
one  (m  =  1)  ,  meridional  index  one  (n  =  1)  is  forced.  The  results 
from  these  cases  are  fairly  similar  so  only  the  longest  wave 
case  (case  1)  is  analyzed  in  detail. 

Two  classes  of  cases  are  examined  using  the  steady  state 
values.  In  the  first,  the  magnitude  of  the  heating  function 
is  changed  and  the  model  is  integrated  from  the  steady  state 
values.  In  the  second,  the  nagnitude  of  the  steady  state 
values  are  changed  and  then  used  as  initial  conditions  for 
further  model  integrations.  In  the  these  cases  the  heating 
function  is  the  same  as  the  one  originally  used  to  bring  the 
model  to  a  steady  state. 
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IV.  EESDLTS 

Due  to  the  similarity  of  the  six  linear  cases  only  cases 
1,  5  and  6  will  be  discussed.  These  three  cases  represent 
the  three  zonal  wavenumbers ( 1 ,  2  and  3)  and  the  tvo  basic 
meridional  structures  possible  in  the  model.  Cases  1  and  6 
are  symmetric  in  their  meridional  structure  while  case  5  is 
anti-symmetric.  The  primary  tool  used  to  analyze  these 
cases  was  the  harmonic  dial.  These  diagrams  are  polar  plots 
of  phase  and  amplitude  for  a  given  spherical  harmonic.  This 
plot  has  the  advantage  that  the  quasi-sta tionary  and  tran- 
sient nature  of  the  waves  are  easily  seen.  Fig.  4.1  is  a 
good  example  upon  which  the  merits  of  the  harmonic  dial  can 
be  explained.  In  this  figure  the  phase  and  amplitude  for 
the  temperature  wave  for  case  1  are  plotted  every  12  h.  The 
elapsed  time  since  the  start  of  the  integration  is  indicated 
every  5  days.  For  an  ideal  case  where  there  is  a  stationary 
and  a  transient  component  one  would  expect  to  see  a  circular 
pattern  which  is  confined  to  a  particular  quadrant.  This  is 
because  as  the  transient  wave  moves  through  the  stationary 
wave  it  will  come  in  and  out  cf  phase  with  the  stationary 
wave.  The  maximum  amplitude  occurs  when  the  transient  and 
stationary  waves  are  completely  in  phase  and  the  minimum 
amplitude  occurs  when  the  waves  are  180°  out  of  phase.  A 
rough  estimate  of  the  phase  speed  of  the  transient  wave  can 
be  obtained  by  observing  the  change  in  phase  with  time  as 
plotted  on  the  diagrams.  The  direction  of  propagation  can 
be  determined  by  observing  the  direction  of  the  phase  change 
on  the  dial.  A  counterclockwise  change  in  phase  indicates 
eastward  propagation,  a  clockwise  change  indicates  westward 
propagation.  In  Fig.  4.  1  the  phase  speed  of  the  temperature 
wave  at  a  =  .925  can  be  easily  estimated.  The  period   of  the 
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wave    is    21-5  days,       thus    the      phase    speed    in   degrees    per    day 
is   36  0/21.5=16.74.  The    direction    of    propagation      is    east- 


PHASE  VS  AMPLITUDE 
FOR  TEMPERATURE 


Tn 


<J  —  .925 


Figure   4.1        Harmonic  Dial  for    Case    1    Temperature   Component. 

ward.  Another  important  point  to  note  here  is  that,  given 
the  uniformity  of  the  change  in  phase,  the  phase  speed  of 
the  transient  component  is  constant.  The  harmonic  dials  for 
surface   pressure,         temperature   and    the    zonal(U)  and    meri- 

dional (V)  wind  components  for  the  lowest  model  level  are 
shown  in  Figs.  4.2  -  4.7  .  Only  those  components  which  have 
significant   amplitude  are      shown   in    these    figures.  In    all 

cases  the  pressure  and  winds  waves  show  a  much  more  irreg- 
ular pattern  than  the  temperature  wave.  These  irregulari- 
ties might  be  due  to  the  presence  of  a  number  of  transient 
modes  components  with  periods  much  less  than  that  indicated 
by  the  overall  pattern.  Table  II  gives  the  estimated  period 
(in  days)  and  phase  speed  (in  °/days)  for  the  temperature, 
surface    pressure    and      the   V    and   U      wind      component    waves    for 
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Figure   4.2        Harmonic   Dials    for   Temperature 
and  Surface  Pressure   for   Case    1. 
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Figure   4.3         Harmonic   Dials    for   Meridional  (0) 
and    Zonal (D)     Wind   Case    1. 
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Figure  4.4    Harmonic  Dials  for  Temperature 
and  Surface  Pressure  for  Case  5. 
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Figure   4.5         Harmonic    Dials    for   Meridional  (0) 
and    Zonal  (0)     lind   Case    5. 
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Figure  4.6    Same  as  Fig  4.4  for  Case  6, 
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Figure  4,7    Same  as  Pig  4.5  for  Case  6 
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the  three  cases  being  considered.  These  estimates  were 
obtained  from  Figs.  4.1  -  4.7  .  Phase  speeds  were  not 
calculated  for  the  pressure  and  wind  waves  for  case  5  due  to 
the  highly  irregular  nature  of  those  dials.  The  phase  speed 
of  the  temperature,  pressure  and  wind  waves  are  the  same  for 
wave  component  (3,3) .  Such  is  not  the  case  for  the  other 
cases,  though  the  phase  speed  are  similar.  The  computed 
phase  speed  of  the  different  variables  for  case  1  varies 
from  15.6  °/day  for  P  to  18.5  °/day  for  U.  The  other  cases 
have  similar  phase  speeds  for  the  different  components.  The 
general  pattern  which  emerges  from  the  data  is  that  the 
phase  speed  of  the  transient  components  seems  to  be  indepen- 
dent of  wavenumber  and  that  the  waves  are  moving  at  approxi- 

TAB1E  II 
Period (p) ,  and  Phase  Speed  (c)  for  Transient  Components 


Case 

Case 
~"5~ 

Case 

p 

T 
0 
V 

P     c 
22f5   15.6 
21.5   16.7 
19.5   18.4 
21.5   16.7 

P     c 

*  * 

12.5   16.7 

*  * 

*  * 

P     c 
7f5   16.7 

7.5   16.7 

7.5   16.7 

7.5   16.7 

mately  the  speed  of  the  mean  wind (15.5  °/day) .  To  test  this 
hypothesis  additional  integrations  were  accomplished  holding 
all  other  factors  the  same  but  increasing  the  speed  of  the 
mean  wind.  The  results  of  these  integrations  showed  an 
increase  in  phase  speed  of  the  component  waves  corresponding 
to  the  increase  in  the  mean  wind  speed 

A.   STEADY  STATE  CASE 

Since  the  purpose  of  this  study  was  to  see  how  sensitive 
the  steady  state  solutions  are  to  errors  in  the  forcing  and 
initial  conditions   it  is   important  that   the  steady   state 
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solutions  of  the  model  be  examined  to  insure  that  they  are 
indeed  steady  and  that  these  sclutions  are  consistent  with 
known  theory. 

To  ottain  steady  state  solutions  the  basic  equations 
were  integrated  up  to  30  days.  It  is  obvious  from  Pigs. 
4.1-4.7  that  after  a  rather  lengthy  integration  the  model 
atmosphere  has  yet  not  reached  a  steady  state  and  that  in 
fact  the  transient  modes  identified  in  the  last  section  are 
clearly  evident  throughout  the  period  of  integration. 
Though  there  are  indications  that  the  solutions  are  slowly 
converging  as  indicated  by  the  begining  of  spiral  patterns 
in  Figs.  4.1-4.7  a  final  steady  state  would  not  be  reached 
for  some  time.  However,  by  averaging  over  the  period  of  the 
transient  wave,  it's  effect  can  be  averaged  out.  The 
remaining  values  are  the  steadj  state  values.  The  above 
procedure  was  performed  on  case  1.  To  insure  that  the  part 
of  the  solution  which  remained  after  averaging  was  truly  the 
steady  state  solution  another  integration  was  performed 
using  the  averaged  values  as  the  initial  conditions. 
Harmonic  dials  were  again  constructed  for  each  of  the  vari- 
able wave  components.  If  the  fields  obtained  by  averaging 
were  truly  the  steady  state  solutions,  one  would  expect  to 
see  no  phase  or  amplitude  change  with  time,  i.e.  the 
harmonic  dial  would  collapse  to  a  single  point.  Figs.  4. 8 
and  4.9  contain  the  harmonic  dials  from  a  20-day  integration 
initialized  with  the  averaged  values  for  case  1.  The 
temperature  wave  is  almost  completely  stationary,  but  the 
pressure  and  velocity  components  have  small  growth  and  move- 
ment. The  fact  that  the  pressure  and  wind  component  waves 
are  not  completely  stationary  is  another  indication  that 
there  are  transient  modes  in  these  waves  which  are  not 
contained  in  the  temperature  wave.  Consequently,  the  aver- 
aging which  was  performed  over  the  period  of  the  temperature 
wave  did  not  average  out  all  of  the  transient  parts  of  the 
pressure  and  wind  waves. 
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Figure   4.8        Same    as    Fig    4.2   for    Steady 
State  Initial    Conditions. 
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Figure  4.9    Same  as  Pig  4.3  for  Steady 
State  Initial  Conditions. 
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It  appears  from  the  above  integrations  that  solutions 
obtained  by  averaging  over  the  period  of  the  transient 
component  are  very  close  to  steady  state  values.  In  order 
to  better  understand  these  solutions  they  will  be  compared 
to  the  steady  state  solutions  from  a  simplified  analytic 
model.   This  model  will  be  derived  in  the  next  section. 

B.   THE  ANALYTIC  HODEI 

The  predictive  system  of  eguations  for  planetary  scale 
motion  in  Z=-ln  (P/P)  coordinates  are: 

D+||  -Z  =  0  (4.1) 

|-(|±V^v|£^rCZ>+Zr»  =  Q(Z)  C-2) 

V  =  ikx^  ^-3> 

where 

D  =  v^V 

Z  =  -p/p 

4>  =  $(Z)+<f>'  (x,  y,  Z) 

3^  3Z  T  vCp    H  3Z; 

H2  =  RT/g 

r.'"=  HlS(S-  +  i  IT) 

T   Cp    H  3Z    ■ 

Eg.  4.1  is  the  continuity  eguation,  Eg.  4.2  is  the  first  law 
of  thermodynamics  and  Eg.  4.3  is  the  geostrophic  relation. 
These  eguations  were  obtained  ly  scaling  the  general  taro- 
clinic  eguations  (Haltiner  and  Williams,  1980)  with: 
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l/a~1      R      <.01    e  -  100      fla/f  ~1 

Where   L    is   the    horizonal    scale,      R      is    the    Rossby    Number  and  e 

is    a      rotational  Proude      number.        Assuming      a    steady      state 


3_ 
droping    therr(Z)    term,    because    it    is    small,    Eg.    4.2    becomes: 


atmosphere     ( £_  =  q     )       in    solid    rotation       (U=  w a    cos  V     )       and 


The    meridional    and    zonal   components      of    the    geostrophic   wind 
in   spherical  coordinates   are: 


3<j>'  (4.5) 


2aI]cosYsinT    d\ 


U    =    - 


8(f)'  (4.6) 


2a2flsin^    9T 


Using   Egs.      4.5   and    4.6      the   divergence    in    spherical   coordi- 
nates  can   be  written: 


n    _    20, cos^  BV  (4.7) 

2 aft  sin*  f 

Now    Eg.    4.1    can   be    written: 

ii  _7_      ,      1  ii'        n  (4-3) 

3Z      *      a^'sinzy    3A 

Assuming    solutions   of  the   form: 

V     =    $(Z)eimXPm(n  (4-9) 

n 

Z    =    W(Z)eimApm(H')  ^-1°) 

n 

Q    =    N(Z)eimAP™(Y)  <4-1l> 

Eg.    4.4    and    4.8    become: 

a>im$(Z)+rW(Z)     =    N(Z)  (4-12) 
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Eg.       4.12   and   4.13    can   be      combined    to    obtain   a    single    egua 
tion    for   W 


a)a2^sin2Y(W"-W  )+TW    =    N  (4«14) 

Eg.  4.14  is  just  a  2nd  order  ncnhomogenous  ordinary  differ- 
ential   eguation    for    "H   whose    general    solution   is    of    the    form: 

w    =    C1W1+C0W0+W 
11       11       p 

Where  F]_  ,  >f?  are  the  solutions  to  the  homogeneous  part  of 
Eg.  4.14  and  W p  is  the  particular  solution.  Assuming  N  (Z)  = 
N0e"a  ?   the  solution    to  Eg.    4.14   is 


-    Cler  +  Z  +  C2er-Z+Ae-a2Z    =    eZ/2  (  C^  ^Pe"^  Z)  l«-15> 


W    =    CleI,  +  i  +  C2er-^Ae-u    "    =    ^"(C^2 

where 

r±    =    y   i/lM-r/coazfisiny      =    |   +  q2 

2 
a  =   seme  positive    number 

A    =    N0/(a)a2sin2nalt+a2  )+D 

There  are  two  possible  solution  cases  for  the  homogeneous 
part  of  Eg.  4.14.  If  g2  is  positive  the  solution  consists 
of      an      exponentially      growing    and      decaying      part.  After 

applying  the  boundary  conditions  w(0)=0  and  finite  energy  at 
Z=°o,    i.e.  (  p|  <J>I  2    )    =    e-z|<H2      the  solution    becomes 

W(Z)    =    A(e-a2Z-eZ/2e-<l2Z)  ^  <"-16> 

An  expression  for  $ (z)  can  be  obtained  by  substituting  Eg. 
4. 16    into   Eg.    4.  13 
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2 


$(Z)    =    -a2fisin24'A{erZ(r-l)+e"a    Z(l+    ct2)}  (4.17) 

If  g2  is  negative  the  solution  will  be  oscillatory  in  Z  and 
of    the    form 

W(Z)    =    eZ/2CC1e-illZ+C2feillZ)+Ae-a2z  <4'18> 

where 

u    =    /r/ooa2fisin24/-    1/4  (4.19) 

After  applying  the  radiation  boundary  condition  at  Z=°°i.e. 
w (Z)  6  (Z) >0  and  the  lower  boundary  condition  w  =  0  at  Z  =  0  Eg. 
4. 18    becomes 


W(Z)    =    A{(|i-y)e(1"2+i^)Z-i(-l-a2)e-a    Z}  <"-20) 

The    corresponding    expression    for         (z)     is 

*(Z)    -    a2fisin2yA{(|i-u)e(^2  +  ili)Z-i(-l-cc2)e-a2Z}  {U-21) 

Egs.  4. 16,4.17  and  Egs.  4.20,4.21  can  be  substituted  into 
Egs.  4.9  and  4.10  to  obtain  expressions  for  <J>and  Z.  After 
taking    the   real    part    the   g2    >    0    solution    are 


d>'     =    a2fisin2W{erZ(r-l)+e"a2Z(l  +  a2)}sin(mnpm  {*m22) 

n 

7    -    a/    -a22      rZ.  ,    ,,Dm  (4.23) 

Z    -    A(e  -e       )cos(mX)P 

n 

The    g2    <    0    solutions   are 

<f>'     =    a2ftsin2YA{eZ/2(|sin(yZ+m\)  (4.24) 

+  Mcos(pZ  +  rnA)+(l  +  a2  )e"a    Zsin(mA)}  Pr 


n 
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Z  =  A{e"a2Z-(eZ/2cos(yZ)  )cos(mA)  (4* 25) 

+(eZ/2sin(yZ)  )sin(mA)  P™ 

If  g2  >  0  the  response  of  the  atmosphere  to  heating  is 
trapped  near  the  surface  and  the  vertical  structure  is  one 
in  which  the  waves  have  no  tilt  with  height.  If  (j2  <  0  the 
perturbation  due  to  the  heating  will  radiate  to  great 
heights  and  the  waves  will  tilt  westward.  In  both  cases  the 
steady  state  represents  a  balance  between  advection,  adia- 
batic  warming (cooling)  and  diatatic  heating  (cooling) .  The 
longitudinal  and  latitudinal  structure,  for  each  case,  will 
be  determined  by  the  forcing  specified.  That  is,  if  waven- 
umber  one  is  specified  in  the  forcing  term,  the  steady  state 
solution  will  have  a  wave  one  form.  In  addition,  since  the 
winds  are  geostrophic  their  structure  will  be  determined  by 
the  structure  of  <J> ' .   Eg-  4.24  can  be  written 

0'  =  -A{(|eZ/2sin(yZ)-eZ/2ycos(yZ)  )cos(mA)  (U*26) 

+(|eZ/2cos(uZ)-yeZ/2sin(UZ)-(l+a2)e"a2Z)sin(mA)}Pm 
which   can   be  written   in   the    form 

(G(Z)cos(nA+H(Z)     )}A(Y)Pm(4M  (4-27) 

n 

This  form  allows  comparison  with  the  phase  and  amplitude  of 
the  model  Legendre  functions.  If  the  model  steady  state 
solutions  really  are  steady  state  one  would  expect  them  to 
at  least  gualitatively  agree  with  the  analytic  solutions 
derived  above.  The  differences  between  the  numerical  model 
atmosphere  and  the  analytic  atmosphere  are: 

1.  A  nondimension  vertical  ccordinate  Z  is  used  in  the 
analytic   model  while   sigma  coordinates   are  used  in 
numerical  model, 

2.  Friction  is  not  included  in  the  analytic  model,  and 
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3.   T(Z)  is  assumed  to  be  constant  in  the  analytic  model 
while  it  is  not  in  the  numerical  model. 
T(Z)  is  not  constant  in  the  numerical  model  because  r(z)   is 
constant.    The  fact  that  T(Z)  is  not  constant  given  r  (z)  is 
can  be  shown  by  examining  the  eguation  for  r (Z) 

r(z)  =^§p+hM>  (*-28) 

After   converting    to    z  coordinates   Eg.    4.28    becomes 

H\g    +    12\  (4.29) 

Since  the  lapse  rate  is  a  constant  6.5  °/Km  in  the  model  the 
second  term  in  Eg.  4.29  is  constant  but  T  decreases  with  z 
so  the  first  term  is  not  constant  thus  r  (Z)  will  not  be 
constant  in  the  model  atmosphere.  However,  since  the 
decrease  in  temperature  is  rather  small  compared  to  the 
mean  temperature  T(Z)  is  nearly  constant  so  that  gualita- 
tive   comparisons  are  still  possible.     In  addition,    the 

vertical  structure  of  the  forcing  in  the  analytic  model  is  a 

— ct2Z 
simple  exponential   form  e      while  in   the  model  it   is  of 

xi.         -e  cr-1  P-Ps  _  .   ,  .  ,  , 

the    form   e~rj  -    e^^    .        However,    given    that        is    chosen   such 

that    the  e-folding    depth    of   e  is      the      same      as    that    for 

P-Ps 
e^TT      then   the      two    functions    will    be   similar   except    for    the 

upper    levels   where   both    functions      are    small   and    have    little 

effect    on   the    total   solution.         Fig.       4.10a    is    a    plot    of    the 

numerical      model      forcing      vertical    structure      where      P        is 

3  s 

assumed  to  be  constant  and  egual  to  1000mb.  This  is  not  a 
bad  assumption  because  the  model  contains  no  topography  and 
the  maximum  pressure  perturbations  are  only  a  few.  millibars. 
This  assumption  also  allows  gualitative  comparisons  between 
variables  in  and  P  coordinates.  The  e-folding  depth  here 
occurs  at  P=700mb.  For  P=700  ,  Z=-ln  (700/1 000)  =-ln  (.  7)  . 
Thus,    to   have   the   same   folding    depth      az       must    be    choosen    to 
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Figure  4.10   Heating  Vertical  Structure 
(a)  Numerical  model  and  (b) Analytic. 
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be    equal   to    1/(-ln.7).         Fig.       4.  10b      is      a      plot      of      e-0*2 
converted   to  P      coordinates.  Note    that    the      plots    are    very 

similar  up  to  400mb.  Figs.  4.11a  and  4.11b  are  plots  of 
G(Z)  and  H  (Z)  vs  pressure  usirg  the  following  estimate  for 
the  various  parameters  in  Eg.  4.19:  w  =(20  m/s)/a,  r~  (Z) 
H2N2,  H=10km/  N2=10-*s-i,  a=6.37  x  106  m  and  ^  =  7.292  x 
10-ss-i  .  Figs.  4.12  and  4.13  contain  plots  of  the  steady 
state  model  solutions  for  the  neridional  (V)  and  zonal  (U) 
wind    vs   sigma    for   case    1.  Since    the    winds   are    geostrophic 

the  vertical  structure  U  and  "V  can  be  compared  with  the 
vertical   structure      of    <}>.  Note   the      striking    similarities 

between  the  two  sets  of  figures.  Both  amplitude  plots  show 
a  slight  decrease  in  the  lower  layers  and  an  exponential 
increase  near  the  top  of  the  atnosphere.  Also,  given  Figs. 
4.14,  4.15  and  4.16,  it  is  obvious  that  the  model's  steady 
state  solutions  have  the  longitudinal  wave  one  form  speci- 
fied in  the  case  1  forcing  term.  These  figures  are  hori- 
zontal cross-section  for  vorticity,  meridional  velocity  and 
temperature  near  the  top  and  bottom  of  the  model  atmosphere. 
Contour  intervals  are  indicated  at  the  top  of  each  cross- 
section.  The  vertical  scale  at  the  bottom  of  the  figures  is 
used  to  represent  topography,  which  is  not  included  in  this 
model.  It  is  clear  that  the  steady  model  solutions  are 
qualitatively  consistent  with  the  analytic  solutions  and 
that  the  method  used  to  determine  these  values  was  a  sound 
one.  The  guestion  might  now  be  asked  why  are  the  model 
solutions  oscillatory  in  nature?  Is  it  not  possible  that 
they  could  be  of  the  exponential  decaying  type?  These  ques- 
tions can  be  answered  by  examining  Eg.  4.19  more  closely. 
Given   that 

r/a^sin2^    >   i  (y-30) 

4 
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Figure  4.  12    Model  Vertical  Amplitude 
and  Phase  Structure  for  Case  1. 
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the  solution  will  be  oscillatory.  In  the  model,  oi  is  speci- 
fied so  that  the  zonal  wind  is  equal  to  20  m/s  at  the 
equator  so  that  U=  co  a  cos  =  U=coa  at  =  0  thus  co  =  U/a.  Given 
r~H2  N2#  H~10  km,  N2  =10"«s-i,  a=6.37  x  10*  m,  Q  =7.292  x 
10-5s-1  one  can  try  to  solve  fcr  a  latitude  where  the  solu- 
tion will  no  longer  be  oscillatory.  After  using  the  above 
values    Eg.    4.30    becomes 

sinT  =20.84  (*-31) 

which  has  no  solution.  Thus  there  is  no  point  where  the 
solution  would  not  be  oscillatory.  To  obtain  non- 
oscillatory  solutions  one  would  have  to  increase  the  zonal 
wind  speed  to  almost  100  m/s.  Using  these  steady  state 
solutions  it  is  now  possible  to  examine  the  sensitivity  of 
these  solutions  to  changes  in  the  forcing  and  initial 
conditions. 

C.   SENSITIVITY  EXPEBIUENTS 

Using  the  steady  state  solutions  from  case  1  as  initial 
conditions  two  additional  experiments  are  carried  out  to 
determine  the  sensitivity  of  the  planetary  waves  to  'errors' 
in  the  forcing  or  initial  conditions.  In  the  first  experi- 
ment errors  are  introduced  in  the  forcing  field.  In  fact  a 
100%  error  was  introduced  by  turning  off  the  forcing  alto- 
gether. The  model  is  then  integrated  out  to  20  days.  Figs. 
4.17  and  4.18  contain  the  harmonic  dials  for  the  above  inte- 
gration. Note  that  for  the  most  part  the  once  quasi- 
stationary  solutions  are  now  almost  purely  transient  since 
the  phase  through  the  entire  range  of  phase.  In  the  previ- 
ously observed  pattern  the  stationary  plus  transient  nature 
of  the  waves  was  indicated  by  a  change  of  phase  and  ampli- 
tude which   occured  in   a  particular   quadrant  of   the  dial. 
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The  estimated  phase  speeds  obtained  from  Figs  4.17  and  4.18 
show  these  transient  waves  to  again  be  moving  at  the  speed 
of  the  mean  wind.  The  fact  that  the  waves  are  pure  tran- 
sient now  is  not  surprising  since  the  heating  is  no  longer 
present  to  fix  the  wave  to  a  particular  location. 

In  the  second  experiment  errors  are  introduced  into  the 
rotational  part  of  the  steady  state  wind  by  reducing  those 
components  by  20%.  The  model  is  again  integrated  out  to  20 
days.  Figs.  4.19  and  4.20  contain  the  harmonic  dials  for 
this  integration.  From  these  figures  one  can  see  that  the 
only  wave  which  is  moving  eastward  is  the  temperature  wave. 
All  the  other  component  waves  are  moving  westward.  The 
temperature  wave  has  a  phase  speed  of  about  12  °/day  at  the 
equator  which  is  about  25%  slower  than  the  mean  wind.  It  is 
evident  that  there  are  a  number  of  high  frequency  components 
present  in  the  other  component  waves  which  are  not  present 
in  the  temperature  wave.  Since  there  is  no  simple  circular 
pattern  for  the  wind  and  pressure  waves  it  is  difficult  to 
obtain  the  period  of  the  transient  part  of  the  solution. 
The  exact  phase  speed  of  these  wave  is  not  important.  The 
important  thing  to  note  here  is  that  the  waves  are  moving 
slowly  westward. 
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Figure  4.17    Same  as  Fig  4.2  fcr  the  Case  with  No  Heating. 
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Figure  4.18    Same  as  Fig  4.3  fcr  the  Case  with  No  Heating. 
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Figure  4.19    Same  as  Pig  4.2  for  the  Case 
with  Rotational  Wind  Errors. 
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Figure  4.20    Same  as  Pig  4.3  for  the  Case 
with  Eotational  Wind  Errors. 
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V.  DISCUSS IOH  ANE  CONCLUSIONS 

The  results  of  the  previous  section  are  best  discussed 
in  terms  of  normal  mode  theory.  A  number  of  investigators 
(Geisler  and  Dickinson,  1976,  Daley  et  al,  198  1,  and 
Kasahara  1976)  have  used  this  method  to  examine  the  tran- 
sient part  of  planetary  waves.  The  method  basically 
consists  of  solving  the  linearized  set  of  equations 
describing  oscillations  of  a  stratified  resting  atmosphere 
on  a  spherical  earth  by  separating  the  equations  into  a 
latitudinal  structure  equation  and  vertical  structure  equa- 
tion. The  vertical  structure  equation  determines  the 
vertical  modes  of  oscillation.  The  separation 
constants (eigenvalues)  from  the  vertical  structure  equation 
can  be  used  to  find  the  eigenfunctions  of  the  latitudinal 
structure  equation.  The  latitudinal  structure  equation  for 
a  given  eigenvalue  is  identical  to  a  fundamental  equation  of 
a  free  surface,  if  the  eigenvalue  is  replaced  by  the  uniform 
depth  of  the  fluid.  This  equation  was  first  obtained  by 
Laplace  in  his  study  of  free  oscillations  of  shallow  water 
over  a  rotating  sphere.  For  this  reason  the  horizontal 
structure  equation  is  often  referred  to  as  Laplace's  tidal 
equation  and  the  eigenvalues  are  referred  to  as  equivalent 
depths.  The  mode  associated  with  the  largest  equivalent 
depth  is  referred  to  as  the  external  or  barotropic  mode. 
This  is  because  this  mode  has  a  vertical  profile  of  hori- 
zontal divergence  with  the  same  sign  throughout  the  atmos- 
phere. The  other  modes  associated  with  smaller  equivalent 
depths  are  referred  to  as  internal  modes.  These  modes  are 
oscillatory  in  the  vertical.  Por  a  given  equivalent  depth, 
zonal  wavenumber  and  meridional  mode  there  exist  three  sets 
of  eigenf requencies  and  eigenfunctions.     Two  are  eastward- 
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and  westward-propagating  gravity  waves  and  the  third  is  a 
westward  propagating  Rossby-Haurwitz  type  wave.  The  phase 
speeds  of  these  waves  are  related  to  equivalent  depth,  meri- 
dional mode  and  zonal  wavenumber.  For  the  Rossby  type  wave 
the  fastest  phase  s  peed  occurs  with  the  largest  equivalent 
depth  and  smallest  meridional  mode  and  zonal  wavenumber. 
The  period  of  this  gravest,  symmetric  zonal  wavenumber  one 
Rossby  mode  is  almost  exactly  5  days.  There  is  observa- 
tional evidence  for  the  existerce  of  this  wave (Madden  and 
Julian,  1972)  and  there  is  also  evidence  that  some  of  the 
planetary  wave  error  in  numerical  models  is  in  part  due  to 
the  spurious  excitation  of  these  large  scale  external  Rossby 
modes  (Daley  et  al,  1981) .  The  following  transient  modes 
which  were  observed  in  the  three  different  parts  of  this 
study  are: 

1.  the  mode  which  moved  at  the  speed  of  the  mean  wind 

in   the  case   where   the   model  atmosphere   was   being 
brought  to  a  steady  state, 

2.  the  mode  which  moved  at  the  speed  of  the  mean  wind  in 
the  case  where  the  model  heating  was  turned  off, 

3.  the  westward  moving  modes  in  the  case  where  error  was 
introduced  into  the  rotational  part  of  the  wind. 

These  modes  may  be  explainable  in  terms  of  these  external 
and  internal  Rossby  modes.  The  mode  which  is  moving  at  the 
speed  of  mean  wind  may  be  a  slow  westward  moving  internal 
mode  which  is  overwhelmed  by  the  mean  wind.  The  westward 
moving  modes  in  the  last  experiment  are  most  likely  internal 
modes  as  well  but  are  probably  associated  with  larger  equiv- 
alent depths  than  those  of  the  earlier  experiments.  Since 
phase  speed  for  a  given  wavenumber  and  meridional  mode 
increases  with  increasing  equivalent  depth,  these  waves 
would  have  a  greater  westward  phase  speed  and  would  be  less 
counterbalanced  by  the  mean  wind.  Given  a  large  enough 
westward   phase  speed  one   would   actually  see   a   westward 
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propagation  indicated  on  the  harmonic  dials.  The  low 
frequency  modes  in  the  first  two  parts  of  this  study  were 
most  likely  excited  by  the  discontinuous  nature  of  the 
forcing.  Numerical  studies  by  Geisler  and  Dickinson  (1976) 
indicate  that  even  when  the  forcing  is  gradually  turned  on 
that  low  frequency  Rossby  modes  are  excited.  An  analytic 
study  by  Clark  (1972)  also  showed  that  a  mode  which  moved  at 
the  speed  of  the  mean  wind  would  be  excited  by  a  switch  on 
of  a  vertical  velocity  or  temperature  disturbance.  However, 
this  mode  was  only  cne  of  a  nuiber  of  modes  which  would  be 
excited  and  it  was  also  found  to  decay  with  time.  Thus,  it 
is  possible  that  the  transient  modes  which  were  observed  to 
move  at  the  speed  of  the  meat  wind  are  in  fact  internal 
Rossby  modes.  The  fact  that  internal  Rossby  modes  might 
have  been  excited  when  errors  were  introduced  into  the  rota- 
tional part  of  the  wind  is  consistent  with  the  results  of 
Daley  et  al  (1981)  who  found  that  errors  in  the  rotational 
wind  will  cause  a  spurious  excitation  of  these  modes. 
However,  it  must  be  pointed  out  that  since  the  results  of 
this  study  were  analyzed  using  spherical  harmonics  rather 
than  normal  modes  it  is  impossible  to  confirm  the  above 
hypothesis.  The  question  that  might  now  be  asked  is  why 
wasn't  the  external  iarotropic  node  excited  in  these  experi- 
ments? The  studies  by  Geisler  and  Dickinson  and  Clark  found 
that  this  mode  would  be  excited.  Clark  found  that  this  mode 
would  in  fact  be  the  dominate  mcde.  The  crucial  difference 
between  the  analytic  model  of  Clark  ,  the  numerical  model  of 
Geisler  and  Dickinson  and  the  model  in  this  study  is  that 
this  model  does  not  have  a  vertical  wind  shear  while  the 
aforementioned  models  do.  This  lack  of  vertical  shear  is 
the  most  probable  reason  why  the  barotropic  mode  was  not 
excited.  In  the  absence  of  a  mean  vertical  shear  the 
heating  specified  in  this  model  will  produce  low  pressure  at 
the  surface   and  a   high  pressure   aloft.    The   vorticities 
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associated  with  each  will  be  of  the  same  magnitude  so  that 
the  vertically  averaged  vorticity  vanishes  and  there  can  be 
no  barotropic  mode  response.  When  a  mean  vertical  shear  is 
present  there  is  a  differential  advection  affect  and  the 
upper  and  lower  vorticities  dc  not  exactly  cancel,  thus 
allowing  the  possiblity  of  a  barotropic  mode 

Given  the  slow  phase  speed  of  the  transient  modes 
excited  in  this  study  by  'errors'  in  the  forcing  and  initial 
conditions  it  does  not  appear  that  planetary  waves  are  very 
sensitive  to  errors.  However,  this  apparent  lack  of  sensi- 
tivity is  most  likely  due  to  the  simplified  conditions  used 
in  the  model  rather  than  to  the  inherent  nature  of  planetary 
waves.  Specifically,  the  lack  of  wind  shear  prohibits  the 
excitation  of  the  fast  external  Rossby  wave.  The  external 
mode  has  a  fast  phase  speed  (  72  °/day)  and  if  excited  it 
could  rapidly  propagate  throughout  the  region  of  integra- 
tion. 

The  results  obtained  by  heating  in  this  study  might  also 
have  been  obtained  by  including  topography  along  with  a 
stronger  surface  friction  term.  Future  studies  should 
investigate  the  sensitivity  of  planetary  wave  to  errors  in 
topography  as  well  as  heating.  In  addition,  future  studies 
will  have  to  use  a  more  realistic  vertical  wind  profile. 
Due  to  the  lack  of  a  realistic  wind  profile,  which  prohibits 
the  excitations  of  the  external  barotropic  mode,  any  conclu- 
sions about  the  sensitivity  of  planetary  waves  to  initial 
conditions  and  forcing  must  be  delayed  until  more  comprehen- 
sive studies  are  completed. 
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